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ABSTRACT 


A general purpose software package was developed to 
perform nonlinear constrained optimization of user-defined 
engineering design problems of significant complexity using 
desktop computers. The package, designated Microcomputer- 
based Design Optimization Tool (MDOT), will accept nonlinear 
functions of up to ten variables, which may be bounded, with 
as many as fifty constraints. It was implemented on a 
Hewlett-Packard Model 85 microcomputer with 32 Kbytes of 
random access memory. 

MDOT employs the method of feasible directions for con- 
Sunatmed Optimization, and a variable metric method for 
Umeonstrarined functions. [t 18 interactive, provides for 
monitoring the optimization progress, and can be interrupted 
to restart from a new point in the design space. Typical 
applications of MDOT are in the design of machine com- 


ponents, composite laminates, and piping systems. 
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NOMENCLATURE 


The nomenclature defined here is that used in the text. 


Definitions of parameters associated with the program code 


accompany the code listings in Appendix B. Boldface charac- 


ters denote vectors or matrices. 


a 
D 
2 
G 
H 
1 


m 


n 


p 


q 


One-dimensional search step length 


Inverse Hessian approximation update matrix 


Objective function 


Vector 


Ofer Mequality constraints 


Hessian matrix, or an approximation to its inverse 


Vector 
Number 
Number 


Vector 


of design variable lower bounds 
Otsrmeduality —constraines 
of design variables 


used in constructing D 


Iteration number 


search 
Sealer 
Sealer 
Vector 
Sealer 
Vector 


Vector 


direction vector 

used in constructing D 

used in constructing D 

of design variable upper bounds 
used in constructing D 

of design variables 


used in constructing D 


objective function gradient vector 





ACKNOWLEDGEMENT 


The author gratefully acknowledges the technical assist- 
ance rendered by Professor G. N. Vanderplaats, the equipment 
support provided by Pat Carroll and Vee Masuero of Hewlett- 


Packard, and the patience of his lovely wife, Daphne. 





I. INTRODUCTION 


Rea OBJECTIVE 

This thesis presents, and describes the development of, 
a computer software package: "Microcomputer-based Design 
Mecimization Tool" (MDOT). The motivation for this work 
stemmed from the lack of available general purpose programs 
capable of performing nonlinear constrained optimization of 
engineering design problems of significant complexity using 
desktop computers. 

In a more general sense, MDOT is intended to help focus 
attention on the versatility and computational power of 
microcomputers. These machines are a potentially very 
valuable resource which is just beginning to be tapped by 
the engineering design community. 

The remainder of Chapter I is devoted to an overview of 
where microcomputers stand in engineering design , where 
MDOT stands amid the optimization software currently avail- 
able, and the implementation of MDOT. In Chapter II, a 
general description of optimization concepts and methods, 
and their application to engineering deSign, is presented. 
In Chapter III, the program development of MDOT is de- 
scribed, and flowcharts of the algorithms coded are pro- 
vided. In Chapter IV, the test problems which were used to 


validate MDOT are described, along with the solutions 





obtained. Chapter V is a brief summary. Appendix A is the 
MDOT user manual. Appendix B contains an annotated listing 


of the MDOT program code. 


B. OVERVIEW 

There are desktop computers available today with memory 
size and computational speed in excess of those of the 
mainframes of just a few years ago. The fact that their 
capabilities are not yet being fully exploited in the day- 
to-day process of engineering design can be attributed in 
part to the lack of available software. As Falk [Ref. 1: 
p.50] observes, "...engineers...-have little time or pa- 
tience to do computer programming." Even among those engin- 
eers who have the time and patience, there persists a reluc- 
tance to program on microcomputers because of a perceived 
lack of general purpose utility or under-estimation of the 
capability of these machines. 

Desmen Optimization 1S a coneept which, similar to the 
desktop computer, has received "mixed reviews" from the 
engineering disciplines. While there are few who would ques- 
tion the virtue of seeking the "best" solution to a problen, 
there are many who are reluctant to relinquish to a computer 
what they see as the engineers! proprietary decisions in the 
design process . Our whirlwind courtship of computer aided 
design (CAD) is being tempered somewhat by a counter trend 


back toward "human aided design". 
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It was within the framework of these two ideas that the 
development of MDOT was undertaken; not only to make avail- 
able a useful interactive design optimization program, but 
to demonstrate that a powerful general purpose problem sol- 
ver can be implemented in a microcomputer. Little knowledge 
of programming is required of the MDOT user. Problem entry 
and program execution are convenient. The interactive fea- 
tures of the code permit the design engineer to keep in 
close touch with the progress of the problem solution and to 
interrupt program execution to make parameter adjustments 
based on engineering judgement. 

The chief advantage that microcomputers enjoy over main- 
frames is their low cost. Small computers are typically 
Dumenased Outright, so that their use incurs no additional 
expense for connection or run time. Their major disadvantage 
is comparatively slow computational speed, but it is doubt- 
ful that engineering design ever progresses so rapidly as to 
make it imperative that a solution be obtained in seconds 
rather than minutes. In any case, this disparity is rapidly 
disappearing with the development of 16 and 32-bit micro- 


processor-based desktop computers [Ref. 2: p.2]. 


C. OPTIMIZATION SOFTWARE CURRENTLY AVAILABLE 
There are several powerful general purpose optimization 
programs available, such as COPES/CONMIN [Ref. 3], which can 


deal with a wide range of design problems. These programs 





must reside in a mainframe computer, and their use can be 
cumbersome, especially for the occasional user. At the other 
extreme are those codes developed for use in computers with 
limited memory. Typically these are special purpose programs 
employing zero order or simple first order methods, such as 
random search or steepest descent, capable of handling only 
relatively small problems. They are convenient, but of lim- 
ited usefulness. 

The gap between these two categories requires that opti- 
mization of the great number of general design and analysis 
problems which are ona scale that could easily be handled 
by small computers be done on a mainframe or not at all. 
This often leads to overmodeling, wherein a relatively 
Simple problem is unnecessarily made more complicated in 
order co more fully utilize the machine capability or to 
justify the expense of computer services. MDOT was developed 
specifically to bridge this gap. 

MDOT provides the design engineer with a convenient tool 
for optimization of nonlinear problems in up to ten bounded 
independent variables subject to as many as fifty inequality 
@enstraints. All that is required is Bee to a desktop 
computer, and today there are certainly few engineers who 


lack this. 





D. IMPLEMENTATION OF MDOT 

The program development for MDOT was done on a Hewlett- 
Packard model 85A microcomputer, which is built around an 8- 
bit microprocessor. This particular machine is so often used 
as a data aquisition system controller that its stand-alone 
computational capability may frequently be overlooked. A 
versatile, engineering oriented computer with high machine 
precision, it is nonetheless on the low end of the memory 
size scale with just 16 Kbytes. The computer's capabilities 
were enhanced by the addition of a 16 Kbyte memory extension 
module and three read-only-memory modules: the matrix, ad- 
vanced programming, and printer/plotter ROMs. As configured, 
there were just over 30 Kbytes of memory available for 
programming. 

MDOT was written in HPBASIC, which differs in some re- 
spects from standard BASIC. As such, in its present form 
MDOT is limited to use in the HP series 80 computers. With- 
Sucemmenmeditriculty, though, the code could be translated 
and run on almost any available hardware. Additional com- 
ments concerning transferrability of the code are presented 


i eemapuer LIT. 
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Eis OPTIMIZATION 


A. APPLICATION TO ENGINEERING DESIGN 

Design has been described as the creative process 
through which the engineering profession develops devices, 
processes and systems to fill the needs of man [Ref. 4: 
p.-170]. A "need" must first be defined in terms of specific 
requirements which the design must meet. Then the engineer, 
drawing on available resources, synthesizes proposed solu- 
tions to meet these requirements. Many such designs may, and 
usually do, exist. Thus there arises the subproblem of 
finding the best of these designs, and the inherently iter- 
ative nature of design. 

imaegirclonally, the solution to this subproblem was 
sought through comparitive analysis of a reasonable number 
of alternative designs; a tedious and expensive procedure 
for problems of even moderate complexity. The recent devel- 
opment of a broad range of very useful CAD software has made 
it possible to remove a good deal of the tedium, and, per- 
haps to a lesser degree, the expense of engineering design. 
For the most part, though, these tools have made no fun- 
damental change in the approach taken to solve the design 
problem. What they have done is redefined the phrase "a 


reasonable number of alternative designs". By programming 





the analysis and comparison tasks into a eompucver, the 
engineer is able to consider many more possible solutions in 
the same amount of time. 

Optimization methods are a significant extension of the 
CAD concept in that they enable the engineer to exploit the 
capabilities of the computer over the entire scope of the 
design process. In optimization, the computer is tasked not 
only with analysis and comparison of previously selected 
designs, but with selection of the designs to be considered 
in subsequent iterations as well. Since this intermediate 
design selection can be quite complex, closing the design 
loop in the computer can lead to a considerable savings of 


time and effort in the search for the optimun. 


B. THE NATURE OF THE PROBLEM 

Fundamental to the economical solution of the design 
problem is that it be quantified and formulated mathematic- 
ally to permit conceptual, rather than physical, manipul- 
ation of resources. The design, then, is specified by ass- 
igning values to a set of independent variables which repre- 
sent its physical characteristics. The measure of goodness 
of the design used for comparison with others is expressed 
as some functional relationship between these variables. The 
requirements placed on the design, aS well as the physical 
limitations of the design itself, define a region in the 


multi-dimensional mathematical design space. The design must 
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reside inside this region to be acceptable. In the termin- 
ology of optimization, the measure of goodness is the 
objective function, the requirements and limitations are 
constraints, and designs which fall within the region 
bounded by the constraints are considered feasible. As 
examples, the set of independent design variables might be 
the cross-sectional dimensions of a structural element, the 
Sejective function its weight, and the constraints its 
maximum allowable stress and size limitations. 

In general, the formulation of a design problem leads to 
an objective and a number of constraints, all of which may 
be linear SG Honlinear funtetvions, explicit or implicit in 
many design variables which themselves are subject to limi- 
tations, called bounds or side constraints. Stated mathe- 


matically (Ref. 5: p.9], the design optimization problem is 


to 
Meee 1 Ze : Fx ) objective function 
Bupygect to: G5(X)<0, Se alte inequality constraints 
Wy S49 SUG vcore Side constraints 
Where: Ren X ay Kon ak design variables 
i “Shepley cecil) lower bounds on X 
peur use ued upper bounds on X 


filmer ne: problem formulation it is more convenient to 


define the objective function as a quantity which is to be 


ul 





maximized, such as efficiency or torque, then the above 


statement may simply be modified to read "Minimize: -F(X)". 


C. THE NATURE OF THE SOLUTION 

Optimization is an application of mathematical theory 
concerning identification of the extrema of functions. In 
multivariable calculus, for example, the method of Lagrange 
multipliers is developed, which provides a closed form 
solution for the extremum of a constrained function. While 
useful for demonstrating concepts and developing methods, 
such analytical techniques are not practical for solving any 
but the simplest of problems. Design optimization methods 
involve numerical approximation techniques and iterative 
search schemes. They are ideally suited to, and in fact made 
practical only through, the use of digital computers. 

Many optimization algorithms have been developed around 
Widely varying strategies. Common to most are the three 


basic tasks that make up one iteration of the solution loop: 


1. Selection of a direction in the design space along 
which to search. 

Succ arech for the most improved design in this 
diaeee tion . 

3. Convergence testing to determine when the optimum 


design has been found. 





For unconstrained problems, these tasks are relatively 
straightforward. Addition of a constraint set may, depending 
on the sophistication of the method employed, complicate the 
Pimst two steps considerably. 

Except in the case of zero order methods, selection of a 
search direction involves calculation of partial deriv- 
atives, for which general purpose optimizers use numerical 
techniques, such as finite forward differences. At anv point 
in the design space, the negative of the gradient of the 
objective function indicates the direction in which the 
objective function is most rapidly decreasing. This may not 
be the best direction in which to search, however, if the 
objective function is highly nonlinear or if the design is 
near one or more constraints. Efficient algorithms variously 
employ constraint gradients, Hessian matrix approximations, 
and previous iteration history information in addition to 
objective function gradients to select the search direction. 

Finding the best improved design along a line in the 
specified direction is termed a "one-dimensional search", 
because the objective and constraints are treated as func- 
tions only of the "distance" along this line from the cur- 
rent design point. Techniques employed in the one-dimen- 
Sional search include the golden section and Fibonacci 
methods, polynomial approximations, and combinations of 
these [Refs. 5,6]. If constraints are present, the best 


improved design may not be the point on this line at which 
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pmemongective function is minimized. If a constraint is 
violated, the design is infeasible, so the search algorithm 
must seek the point at which the objective function is 
minimized while remaining inside the feasible region. 

Part of the optimization problem formulation is specifi- 
cation of an initial design point from which to Start the 
solution process. For constrained problems, the possibility 
Sremocs that this initial design will be infeasible. To 
provide for this, the search direction routine must find the 
direction which will yield the shortest path to the feasible 
region, and the one-dimensional search algorithm must allow 
for the possible necessity of increasing the objective func- 
momen order to attain feasibility. 

: Convergence to a global optimum generally cannot be 
prprenteed. Theory provides the Kuhn-Tucker conditions 
necessary for the existence of an optimum, but these are 
neither convenient to evaluate nor sufficient to define 
aeeumality [Ref. 5° pp. 17-20]. In practice, convergence is 
typically considered to be indicated by one or more of the 
following: 
1. Failure to find a search direction which will lead to 
an improved design. 
BeeeGiyven= a direction, failure to find any significant 
search step length to improve the design. 
3. Finding no appreciable design improvement over a 


specified number of iterations. 
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If the possibility of local minima exists, the optimi- 
zation should be repeated from several different initial 
designs. For constrained problems, the optimization process 
may fail to find any feasible solution, in which case the 
problem must be reformulated. 

An iteration in the optimization solution, then, may be 
summarized. Beginning from the current design point xX, a 
search direction, S49, is determined. Then the one-dimen- 
sional search is conducted to find the "distance", a, along 
S94, which yields the best improved design. The design is 


then updated as 
yA 5 CES 


at which point the objective function is reevaluated and the 
design checked for convergence. 

Optimization algorithmic efficiency and convergence 
behavior are affected by the mathematical characteristics of 
the problem. As numerical methods, they are susceptible to 
ill-conditioning. Truncation and round-off errors, which are 
an unavoidable consequence of the use of digital computers, 
aggravate this. Given an optimizer suitable to the problem 
type, careful problem formulation, as discussed in Appendix 


A, is the best inSurance against poor optimizer performance. 
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Iff. PROGRAM DEVELOPMENT 


ie BASIC CRITERIA 

At the outset of the program development, four basic 
eriteria were established to be met by MDOT: utility, mini- 
mization of required memory, user convenience, and reduction 
of problem run time. At points where conflict existed 
between them, these criteria were prioritized in the order 
listed. Few such compromises were necessary, as the require- 
ments were found to be generally complimentary. 

Mvewuurlity eriterion meant that MDOT should be a 
general purpose optimizer which could be applied to a wide 
range of design problems. Minimization of memory was 
dictated by the limitaions of microcomputers, and affected 
not only algorithm selection and problem size, but many 
aspects of the actual coding as well. User convenience 
considerations drove the development of those portions of 
the code which are interactive, and those involved with 
problem entry and output options. Reduction of problem run 
time was a factor throughout the development, most notably 
iMmenewncorporation of optimization progress display and 


Piewincernupe/restart option. 
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Bae rbOl ALGORITHMS 


1. Algorithm Selection 

Of the many optimization algorithms available, the 
zero order methods, aS well as the simpler of the first 
order methods, were ruled out on the basis of their lack of 
general purpose utility. Others, including linear and quad- 
ratic programming types, were eliminated because of their 
excessive memory requirements [Refs. 7,8]. Finally, the need 
to reduce problem run time while retaining utility lead to 
the selection of two algorithms, each capable of nonlinear 
muULCidimensional optimization. The first is a variable 
metric method which is used in MDOT for unconstrained opti- 
mization, the second is a method of feasible directions, for 
minimizing constrained functions. 

The selection of the one-dimensional Search strategy 
to be employed was driven by the need to reduce problem run 
time. There is a trade-off to be made between the precision 
to which the search step length determination is made and 
the time required for each optimization iteration. Both the 
golden section and Fibonacci search methods can attain very 
precise solutions, but to do so they become computationally 
expensive. In MDOT the one-dimensional search routines were 
designed to seek a leSS precise step length solution in 
order to complete each iteration more quickly. The method 
employed in both the optimizers estimates an initial step 


length based on a reasonable change in objective function 


ee 





magnitude. The golden section ratio is then used to estab- 
lish bounds on the solution, which is finally refined by 


polynomial approximation. 


2. The Unconstrained Optimizer 

From the derivation of the Kuhn-Tucker conditions, 
it is Known that if at some point x* the objective function 
F(X") has a local minimum, then the gradient of the object- 
mye at this point, V F(X”), must vanish and the Hessian 
matrix H must be positive definite. Combined with a second 
order Taylor series expansion of F(X) about some point, say 
X°, near the minimum, these conditions lead to an expression 


for the direction from X° to X* as 


Keo aero) | 9 EGLO) 


In practice, determination of the Hessian matrix by finite 
difference approximation, aS well as inversion of the 
matrix, would be computationally so expenSive as to outweigh 
the theoretical gain in algorithmic efficiency. 

In variable metric methods, information gathered as 
the optimization progresses is used to develop an approxi- 
mation to the inverse of the Hessian matrix, which is then 
used in determining the search direction. As such, these 
first order methods have some convergence characteristics 
comparable to those of second order methods. The algorithm 


for the variable metric method is shown in Fig. 1. 
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The inverse Hessian approximation is initialized as 
ann x n identity matrix. To begin each iteration, the 


search direction is defined as 
So = -H VF(Xd) 


After the one-dimensional search and design update, H is 
modified as 

Hit+l = qd 4 pd 
Where the form of the update matrix D determines which one 
of a family of variable metric methods is being used. In 
general, D is defined as 


Da = StRE ppt + Hol dy(Hdy)? - XCH4ypt+p(Hdy)?) 


Ss \2 S 
where: p= xd - xd-! 
y = VF(X9) -VF(x9-!) 
S = p’y 
a y Hoy 


Two forms of D, and thus two variable metric 
methods, are available in MDOT. The first is the Davidon- 
Fletcher=-Powell method, where wis set equal to zero. The 
second is the Broydon-Fletcher-Goldfarb-Shanno method, with 
w equal to one [Ref. 5: pp. 92,93]. As the convergence 
behavior of a given algorithm can be somewhat problem depen- 
dent, this feature allows the MDOT user to compare the 


results of two variations of unconstrained optimization. 
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The one-dimensional search routine employed by MDOT 
in the variable metric optimizer first finds bounds on the 
unconstrained minimum of the objective function, then 
refines the minimum by a three-point cubic polynomial 


approximation. 


3. The Constrained Optimizer 

The addition of a set of constraints to the optimiz- 
ation problem requires that more sophisticated techniques be 
applied to its solution, particularly in the determination 
of a search direction and the subsequent one-dimensional 
search. As is the case for unconstrained optimizers, itis 
generally the method used to find a search direction which 
distinguishes the different constrained optimization algo- 
rithms. In the method of feasible directions, a search 
direction in which a finite step will reduce the objective 
function is termed useable, while one which will avoid 
constraint violation is called feasible. The direction find- 
ing problem is then formulated as a sub-optimization task to 
determine the best of the possible usSeable-feasible direct- 
ions. MDOT employs the algorithm presented by Vanderplaats 
[Ref. 5: pp. 163-170] for the solution of this sub-problem. 

As shown in the flowchart of Fig. 2, the feasible 
directions optimizer begins aS a Simple Steepest descent 
algorithm, provided the initial design is feasible. Opti- 
mization thus proceeds quickly to a point where one or more 


constraints are encountered. 
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On subsequent iterations the Sub-optimization rou- 
tine is used to determine search directions so as to Satisfy 
the fundamental requirement that the optimum design be 
measrple. The constrained optimizer in MDOT will accept 
initially infeasible designs. In this case, the direction 
and search routines are modified so as to attain feasibility 
as quickly as possible. Thus a feasible direction and step 
length are sought which will overcome the constraint vio- 
lations, even at the expense of increasing the objective 
function. Once inside the feasible region, optimization 
proceeds as before. 

As in the unconstrained case, in the one-dimensional 
search routine employed in MDOT for constrained optimiz- 
peron nw ooundsS on the solution are first established, fol- 
lowed by refinement by polynomial approximation. Here, how- 
ever, the search must be conducted for the zeros of the 
Comsictiatn. functions as well as for the minimum of the 
objective. The step length selected is then the one which 
yields the best feasible design. Provision must also be made 
to ensure the deSign variables remain within their bounds 
(side constraints). In MDOT, if at any time during the one- 
dimensional search a design variable is found to exceed an 
upper or lower bound, it is set equal to the value of the 


violated bound. 
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4. Program Logic 

The relationships between the modules of MDOT are 
depicted in Figs. 3 and 4 for unconstrained and constrained 
optimization, respectively. The main program is named 
"Autost"™ because this signals the HP-85 operating system to 
load and run this program automatically when the computer is 
powered up with the mass storage cartridge inserted in the 
tape drive. All the other modules of MDOT are subprograms 
which are called into memory and executed as needed by 
Autost or another Subprogram. Once entered into main memory, 
a subprogram reSides there for the duration of the optimiz- 
ation unless a "SCRATCHSUB" instruction is executed. 

Following is a brief description of the function of 


each of MDOT's program segments: 


Autost MDOT main calling program 

LOGO Displays introductory (welcome) graphic 
DEF ALT Sets program parameters to default values 
ERO. Problem entry, evaluation of F and G 
CCOnT Gontrot of constrained optimization 

UCONT Control of unconstrained optimization 
ACON Identification of active/violated constraints 
GRAD Calculation of gradients of F and G 
DIRECT Direction finding subproblem solver 
FDSRCH Constrained one-dimensional search 
VMSRCH Unconstrained one-dimensional search 

NEWH Update approximation to Ho! 
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CONV Convergence testing 
PROG Optimization progress information 


TERM Output of results of the optimization 


The PROB subprogram is created by the user by 
editing a skeleton problem entry code which is stored on the 
tape. The edited version is then renamed and stored. Autost 
queries the user for the problem name, which is then common 
wo all subprograms. Autost makes the first call to PROB in 
Srmaer tO select the appropriate optimization control rou- 
tine. Thereafter, PROB is called any time an objective 
eunetton Or Constraint evaluation is required. Both LOGO and 
DEFALT are called by Autost. LOGO generates a simple 
Mwelcome"™ graphics display, and is scratched from memory 
upon execution. DEFALT initializes a number of program para- 
meters to their default values, as defined in Appendix B. 

PROG is called upon completion of each iteration. 
This subprogram generates the user selected optimization 
progress indicators. Options include data and graphics 
fesplays and printed output. Based on the progress inform- 
ation provided, the user may elect to continue the optimiz- 
aieonemerescare MDOI from a different Initial design, or stop 
and reformulate the problem. CONV is also called upon 
completion of each iteration, to determine, based on the 
convergence criteria set by DEFALT, whether the optimum 


design has been found. 
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TERM is called to end the optimization and generates 
the output of results. Termination may be invoked by a 
number of conditions in other than meeting convergence 
criteria. If the optimizer exceeds a specified number of 
iterations, if the components of the search direction vector 
are all essentially zero, if no search step length can be 
found to improve the design, or if there are an excessive 
number of violated constraints, MDOT will terminate. In any 
case, TERM will generate an output message to indicate the 
condition upon which the decision to terminate was based, 
and offer the user the option of editing and restarting the 


program. 


See DAP TATION OF MDOT TO OTHER SYSTEMS 

Sinee MDOT 26 coded in HPBASIC, it is not immediately 
transferrable to hardware other than the Hewlett-Packard 
series 80 desktop computers. Translation of the package, 
either into another version of BASIC or into FORTRAN, is 
certainly a "do-able" project which would significantly 
expand the applicability of MDOT. 

This section highlights those features of HPBASIC used 
in MDOT which would have the greatest impact on this pro- 
ject. They are: variable name assignment, SUBPROGRAMS, 
matrix manipulations, and graphics. 

The limitation of HPBASIC which most decreases the read- 


ability of the code is that variable names are restricted to 
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coueer sone letter or a letter followed by a digit. One 
result of this is that arrays of subscripted variables are 
sometimes used where individual characteristic names might 
otherwise be assigned. A feature of the language which helps 
considerably, not only to overcome this limitation, but in 
programming complicated algorithms, is the SUBPROGRAM. 

Similar to a SUBROUTINE in FORTRAN, the SUBPROGRAM is 
called when needed and variables may be passed either by 
name or by value. In HPBASIC this allows for the use of the 
same variable name to denote different parameters in sep- 
arate program segments. MDOT makes extensive use of SUBPRO- 
GRAMs. This feature is not available in all versions of 
Pisce. saa thout it, the translation of MDOT would be more 
fieebeulen, DUE SELL possible through the use of functions 
ama subroutines, particularly if multi-character Sabie 
names are permitted. 

Matrix manipulation is convenient in HPBASIC. Operations 
pen as Matrix multiplications, transpositions, dot pro- 
ducts, and identifying extreme array elements are all 
accomplished through simple "MAT" statements. This feature 
is available in some of the other versions of BASIC, but not 
in FORTRAN. Without it, additional subprogramming would be 
required to perform these operations. 

Graphics capabilities vary widely from one hardware 
manufacturer to another, as do the coding instructions used 


to execute the displays. it is likely that the graphics pro- 
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grammed into MDOT would require major modification to make 
them transferrable. 

Two more details of the HPBASIC code are worthy of note: 
program flags and the @ symbol. Program flags are built-in 
indicators which can be set to 1 or cleared to 0, and are 
used in MDOT for conditional branching decisions. They could 
easily be replaced by integer variables. The @ symbol is 
used to condense the code and thus conserve memory. It 
Simply separates multiple executable statements on one pro- 
gram line. Without this feature, each statement must have 


its own line number. 


Dee OUuSWITAL FOR FUTURE GROWTH 

Besides its obvious potential for expanded problem size 
if implemented in a computer with a larger memory, there are 
many refinements and additions which could be incorporated 
into MDOT, either to enhance its general purpose utility or 
Bomvailor it to a particular type of problem. Some modifi- 
cations for improved utility might involve coding additional 
algorithms, automatic design variable scaling, and the hand- 
ling of equality constraints. 

MDOT could be customized by modification of the problem 
input subprogram, graphics display, output format, or the 
algorithms themselves. Coupling of MDOT to an external CAD 
or analysis code also presents many interesting and poten- 


tially useful possibilities. One such configuration might be 
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to use MDOT as a Subprogram called to perform optimization 
on localized aspects of a large scale design problem in 


conjunction with a desktop computer CAD system [Ref. 9]. 
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iNew leol CASES 


Validation of an optimization program typically consists 
of testing it on a battery of representative problems to 
Which the solutions are known (Ref. 10]. Based on the 
results of such tests, a number of yardsticks exist by which 
the optimizer is judged relative to others. These can be 
grouped into three categories: stability, robustness, and 
Saueorency (Ref. 6: 9.75]. 

An optimizer is stable if, once a feasible design is 
attained, the objective function remains non-increasing 
until the optimum has been found. A robust optimizer is one 
which elds a valid solution given a poor initial approxi- 
mation. Efficiency refers either to the number of function 
and derivative evaluations required in the solution or to 
the problem run time. These two measures of efficiency are 
closely related if comparing different optimizers run on the 
Same machine, since function and derivative evaluations are 
typically costly operations. 

To these measures of an optimizer's performance, a 
fourth category should be added; that of utility. Given a 
Stable and robust optimizer, there are characteristics in 
addition to its efficiency which should be considered in 


determining its utility. Problem size and type solvable by 
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the program are of fundamental importance, as is memory 
storage required to run it. There are trade-offs involved 
here in program development between the sophistication of 
the algorithms used, speed of convergence, hardware 
capabilities, cost of run time, and user convenience. The 
Mmeinity Of the optimizer is an indication of how these 
trade-offs were made, and involves much more than just 


efficiency. 


A. UNCONSTRAINED TEST PROBLEM 
Among the unconstrained test problems run on MDOT was 


the so-called "banana" function: 
ROH) = UOKy a QOS 2 Oe 2c > ogee 


Mawenhemias atl optimum of F(1,1)=4.0. This Buneeion derives 
its name from the shape of the contours of constant objec- 
tive function (Fig. 5). Although only two-dimensional, the 
banana function is a good test of an unconstrained optimizer 
because the odjective function surface becomes a steep, 
narrow, curved "valley" as the optimum is approached. An 
Met reletens OptLimizer Will tend to "“Zige-zag" in such a 
design space, resulting in slow convergence near the solu- 
tion, while a non-robust optimizer will tend to terminate 


prematurely. 


38 





4.0 


3.0 14 


S, 
OPTIMUM 
oi pg Mes. ? 


oko e =) Unmconstrained Tast Probiaem Design Space 


39 





Results of the performance of MDOT on the banana 


function are summarized below. 


Initial design: Kae =hat 0 
Ol 
KP 2 105 
FO = 10.5 
Optimum: X," = 0.95979 
Xo" = 0.91442 
2 Se One 


MDOT arrived at this solution in about one minute, after 12 


iterations and 66 function evaluations. 


B. CONSTRAINED TEST PROBLEM 

Among the constrained test problems run on MDOT was the 
cantilevered beam problem posed by Vanderplaats [Ref. 3: 
Pmomwemeas illustrated in Fig. 6. The objective function in 
this case is the volume of the beam, for which a theoretical 
optimum of 6603.9 is known. Results of the performance of 


MDOT on this beam design problem are Summarized below. 


Initial design: X,° = 3.5 
clam 
X5 = 50 
FO = 11200 
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Design Variables: b, 4h 
Ob jective Function: Volume 
Canoteaginus: Bending Stress < 20, O00 psi 
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Optimum: X 


T.c201 
X5 = Ve eal ee 


F* = 6637.5 
Gy" = -0.024 
Go" = -0.00656 

% 


G2 = -0.00522 


peel earrived at this solution in about two minutes, after 11 
iterations and 56 function evaluations. 

With refinement of the algorithms, improvement could 
likely be realized in the program performance. All variables 
Perot ware declared "SHORT", which in HPBASIC WGARE they 
are carried to 5 digits. In a machine With just 64 Kbytes of 
memory, this could be changed to "REAL", in which case 9 
digits would be carried, with an attendant improvement in 
the precision of the solution. The number of function eval- 
uations, and thus the problem run time, could be decreased 
by modification of the algorithms such that gradients are 
not calculated in every iteration. Also, the efficiency of 
the one-dimensional Search routines could be improved by 
distinguishing between linear and nonlinear constraint 


funetion's 
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V. SUMMARY 


Optimization is a useful tool in engineering design. The 
desktop computer is the vehicle through which this tool can 
be made Widely available, convenient, and inexpensive. The 
development of MDOT affirms the feasibility of implementing 
a powerful general purpose optimization algorithm in a com- 
puter with limited memory. 

The applicability of MDOT could be expanded through con- 
version to standard BASIC or translation to FORTRAN. It has 
potential for growth in terms of versatility and problem 
pzes and lends itself to tailoring to suit a particular 
class of problem. MDOT could be coupled with a microcomputer 
CAD package to close the design loop in the computer. 

MDOT has been validated by tests on a number of problems, 
both constrained and unconstrained. Its performance is good, 
and could be made better through refinement of the 
algorithms. Specific modifications might involve the one- 
dimenSional search routines and the frequency of gradient 
calculations. 

As microcomputers continue to become more commonplace and 
their capabilities continue to improve, emphasis will shift 
away from the mainframes for the solution of problems which 
are on a scale easily handled by smaller machines. Software 


such as MDOT will both accompany and encourage this shift. 
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AP DENDIEX oA 


MDOT USER MANUAL 


1. INTRODUCTION 

To avoid repetition, references to material presented in 
the preceding chapters are made in this appendix. A useful 
follow-on project would be to assemble a user manual for 
MDOT independent of the background and developmental 
material in the body of the thesis. Details of computer 
operation have not been included here, as it is assumed that 
the user is either familiar with the machine or has access 
to the operating manual. 

MDOT is currently available only on magnetic tape 
cartridge for use in Hewlett-Packard series 80 computers. If 
it is to be implemented in an HP-85A, the machine must be 
configured with four enhancements: a 16 Kbyte memory exten- 
sion module, a matrix ROM, advanced programming ROM, and a 
printer / plotter ROM. No peripheral devices are required, 


nor is extensive programming. 


2. PROBLEM FORMULATION 

Formulation of a well-posed problem, as discussed in 
chapter II, is fundamental to the satisfactory performance 
of an optimization program. First, the design variables must 


be identified. These are the parameters of the problem which 
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the optimizer will be permitted to change in its search for 
the best design. The objective must then be a function of 
these variables, and the minimum of this function is what 
the optimizer will seek. Constraints may be imposed on the 


design in two ways: 


1. Upper and/or lower bounds (side constraints) may be 
specified for any of the design variables. 
2. General inequality constraints may be expressed as 


functions of the design variables. 


Side constraints are explicitly assigned when the initial 
design estimate is entered into the problem subprogranm. 
Inequality constraints must be formulated as quantities 
which are to be less than or equal to zero. Care should be 
exercised to avoid redundant or otherwise unnecessary con- 
straints. In MDOT, an unconstrained problem has neither side 
constraints nor inequality constraints. MDOT has no pro- 
vision for equality constrained problems. 

The objective function can be any characteristic of the 
design expressible mathematically in terms of the design 
variables. It is important to Keep in mind that it is the 
minimum of this function which is sought. If the problem is 
formulated around an objective which is to be maximized, 
then it must be entered in such a way that MDOT will seek to 


minimize the negative of this objective. 
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An important consideration during problem formulation is 
the range of orders of magnitude of the design variables. 
Optimizer performance is best in a design space in which the 
contours of constant objective function are concentric hy- 
perspheres, such that a given change in one variable has the 
same effect on the objective as an equal change in any other 
fmaap le {hef. 11, p-1/]. In practice, this is approximated 
by scaling the deSign variables such that they are all of 
the same order of magnitude, or nearly so. Some optimizers 
G@6ntuis automatically, MDOT does not. 

Selection of the initial design point from which to 
start MDOT will affect its performance and problem run time. 
Any available information which will improve the initial 
approximation should be used. If a constrained problem is 
being entered, a check should be made to enSure the initial 
design falls within the side constraints. Although MDOT is 
equipped to handle initially infeasible designs, convergence 
Will likely be more rapid if the initial design is free of 
violated inequality constraints. In some cases, initial 
Peasipility may be a difficult thing to build into the 
problem formulation. MDOT will display the results of the 
first design evaluation and indicate whether or not it is 
feasible. At that time, the user may elect to proceed with 
the optimization or edit the initial design and restart the 


program. 
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ee) Pp R@meLMPENTRY 
MDOT problem entry is accomplished by editing the PROB 
subprogram. With this module loaded into memory and listed 


on the CRT, modifications are made on the program lines 


noted below. 


Line 10: 
The file name of the subprogram is changed to any 
name up to 6 characters in length, except any of 


those already assigned to MDOT files. 


Eime CC: 
Just after the word DATA, two integers are added, 
separated by a comma. The first is the number of 
design variables, (N1), the second is the number 
Obemmeaguality constraints, (Ne). For unconstrained 


problems N2 is always zero. 


Lines 201-210: 
Initial values of the design variables are en- 
tered, beginning with X1 on line 201 and contin- 
uing, one variable per line. If the problem is 
constrained, each initial value is followed by the 
lower and upper bounds assigned to the correspond- 
ing variable. If no bound is to be specified, the 


biemdeas ftriled by anal". 
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Lines 231-398: 
These lines are available for defining expressions 
to be used in the design evaluation. These state- 
ments will be executed prior to each objective or 
constraint computation. Whenever the design vari- 
ables are used in this, and the remaining sections 
of the subprogram, they are expressed as Xi, 


Pane XO}: 


Lines 400-459: 
These lines are available for defining the object- 
ive function, which must be assigned the variable 


name F. 


Lines 500-9000 
These lines are available for defining the in- 
equality constraint functions, which must be sub- 


scripted variables named G(i), i=1,Ne2. 


Depending on the complexity of the problem, the user may 
elect to use any BASIC programming structures in this sub- 
program. As examples, FOR-NEXT loops, FUNCTIONS, SUB- 
ROUTINES, and even other SUBPROGRAMS could be used in the 
problem formulation. With about 7.3 Kbytes of memory avail- 
able for constrained problem entry, and twice that for 
unconstrained problems, there is space available for con- 
Siderable programming. Variable names used may be any except 


those included in the PROB nomenclature, as defined in 
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Appendix B. If subsecripted variables are to be used, a DIM 
(dimension) statement must be inserted near the beginning of 
the subprogram. When problem entry is complete, the subpro- 
gram is stored under the new file name. 

Listings of the subprograms created for the test cases 
of Chapter IV are included at the end of Appendix B, under 
the names "BANANA" and "BEAM", Comparison of these listings 
to that of PROB, which is the unedited version, will help to 
illustrate the problem entry procedure for both constrained 


and unconstrained cases. 


4, PROGRAM EXECUTION 

MDOT is started either by execution of LOAD "Autost" and 
RUN commands, or by powering up the computer after inserting 
the tape cartridge, as explained in section III.B.4. After 
the "welcome" graphic, the uSer iS queried as to the problem 
mmroenane,. progress display options and output format 
desired. Optimization then proceeds. 

The user may choose to monitor the optimization closely, 
or perhaps not at all. For example, if a constrained problem 
is being run for the first time, the user may want to check 
for rapidly changing design variables so that the option of 
editing and restarting might be exercised. On the other 
hand, if a problem Known to be well-behaved is being re-run 
With relatively minor changes, the uSer may elect to 


Signore’ MDOT until the solution is obtained. 
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Any of the values carried by the DEFALT subprogram, as 
defined in the listing in Appendix B, may be changed. They 
may be permanently modified by editing and re-storing 
DEFALT. Alternatively, a default value may be changed during 
program execution, whenever the edit option is invoked, by 
reassigning its value from the keyboard. If the program is 


restarted, though, DEFALT is recalled. 
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APPENDIX B 


ANNOTATED PROGRAM LISTINGS 


In an effort to conserve memory, MDOT was coded without 
remark statements. The program listings are therefore pre- 
sented in an annotated format to aid in following the logic 
flowpaths. Nomenclature common to two or more modules is 
defined first, then the additional nomenclature unique to 
each module immediately precedes the applicable segment of 
the code listing. Numbers in parentheses in the function 
descriptions refer to line numbers in the associated program 
segment. Parentheses following a variable name in the nomen- 
clature lists indicate vectors, while parentheses enclosing 


a comma indicate two-dimensional arrays. 


COMMON NOMENCLATURE 


A1() Addresses in G() of violated and active constraints 
A2() Values of constraints identified in A1() 

C1c) Vector of integer default values 

eau) Vector of non-integer default values 

D() Current objective function gradient 

D1() Previous iteration objective function gradient 

F Current value of objective function 

FO inp@eralevalue of objective function 
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id] Previous iteration value of objective function 


G() Vector of current constraint values 

Git, ) Gradients of violated and active constraints 
HC, ) Approximation to the inverse of the Hessian 
L() Vector of design variable lower bounds 

N1 Number of design variables 

Ne Number of inequality constraints 

N3 Number of currently violated constraints 

N4 Number of currently active constraints 

P$ Problem subprogram file name character string 
PO-P3 Polynomial approximation coefficients 

Q1 Iteration counter 

Q2 Bunetron evaluation counter 

Q3 Convergence counter 

Q4 Convergence counter 

R 1 Golden section ratio 

R2 Golden section ratio 

a) Search direction vector 

X() Current design 

moc) Initial design 

aC) Previous iteration design 

X9() Working vector of perturbed design variables 
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PROGRAM FLAGS 


Unless otherwise noted, the flags listed are set to a 


value of 1 when the asSociated condition exists. 


FLAG CONDITION 
1 Maximum number of iterations has been exceeded 
2 Number of violated constraints greater than N1i+2 
3 Search vector components are all essentially zero 
4 No appreciable move parameter can be found 
5 Convergence to optimum has occurred 
6 Variable metric algorithm is to be restarted 
7-8 Unassigned 
9 Subprogram is to be exited without execution 
10 Termination has occurred 
14 Progress option select: 1 data display 


2 graphic display 


a8) 





Module: Autost 
galls: Autost calls all other modules of MDOT into 
memory, then scratches those not needed for the 
type of problem (constrained or unconstrained) 
being run. For information transfer, Autost 
Garis: PHOB (edited), UCONT, CCONT, LOGO, and 
DEFALT. 
Function: MDOT main calling program (10-370). 
Interactive problem initiation (1000-1200). 
Program loading (2000-2070). 
Initial design display generation (3000-3130). 
Design/default editing (4000-4260). 
Nomenclature: 
A$ Interactive query response (string) 
S$&A$ Concatenated string to call appropriate optimizer 
D$ Diselay estrone eo ohomeLe™. or 'PNFEASIBLE™) 
B$ tesplay string ("'OC1ret or "C2() 
NO Working variable for editing default values 
N Interactive query response (numeric) 
KO ithe= 0; Iindicates@rmarst problem of run 
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ort FMCG S CT. 2 we 
2® OUrFPION BASE 1 
s@ COM FPtC06] , INTEGER N1,N2,N3,N4,01,92,02,94,01(14) , SHORT 
eee <2) ) 
44 DIM 232049, A$011,0¢010),G403) 
eee ahigim) A9OC10) ,FO,6( 30) L619) ,UC19) , Ne 
4@ INTEGER N,I,Ka 
7@ S$="CONT" 
20 GOSLIE Bawa 
Pomorie “CUI” @ SCRATCHSUE “LoGa" 
14@® CALL "DEFALT" 
119 01,02,03,04,N3,N4,k0=0 
moet l=1 TO 11 
1320 IZTFLAG I 
140 NEXT I 
120 WAIT 1904 @ GEOSuUE 16@0@® 
14@ CLEAR @ CALL Ft ¢ @ ) 
P7o REUQIM XO@C(Ni) 
180 IF N2=0 THEN 219 
V7O Belim is(N2),L(N1),UCN1L) 
2490 GOTL 22o 
=14@ MAT G=ZER(1)@ MAT L=G@ MAT U=L 
Meomermemn e+ € 1,X0()  ,FO,GCdGL¢),UC) >) 
224 GIS sooo 
“40 IF N2=o@ THEN At="LI" ELSE At="" 
PatymealtL Ateoe ( X@C),FO,GC),L4),UC) ) 
~60 CFLAG 10 @ CFLAD 7 
oer © Hite © DISP “Select ontion..." @ ODISF @ DISr * 
1) EDIT/RESTART" | 
Safe € OIsSrF " Zo NESE RUboceMt = € Der Genisrr® ao 
rit” 
“970 INPUT N 
200 ON N GOT S20, 214, 240 
SEeOeskmalrtmnslb Fe C GUT 10H 
324 GBNSLIB 44a 
SeeOeCAct FS ( 1,X8¢),F@,60),L4),U0) ) 
340 GhIstle 3aa@ 
Baw GUT 28 
249 CLEAR 
370 END 


= 
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Lame 
1010 
Lazo 
moat 
1Q3 


strained" 


1040 
LOA 


NEWH" 


1940 


AION" 


197) 


IF F@:@ THEN 1070 
EG 1 
EEEAR © 
Meets ~ @ LTS 
ise he) 
@ INPUT WN 
LIN N GUT 1050, 1040 
SCRATCHSUB “CANT” @ 
@ GOTO 1076 
SICRATCHSUE "CONT" @ SCRATCHSUE 
@ ‘SCRATCHELIE "DIRECT" 
CLEAR @ DISP @ DISe 


Peer bis Woelect prahlem type far this pad 


Loanstrained"” @ DISe @ DIsSpP " —) Lincan 


SC Reem oen Vomit! 2 re RATEH Stina 2 


NEVER iat 


ip 


See eb cwee 


at dems" 


1S 
Layo 
Liga 
mio 
1a, 


1138 
1140 
1150 
1140 


mi / 


L1s@ 
1170 
1200 
~ Oa 
=O1@ 
Ni" 

2O20 
20.30 
ROH" 
O40) 
ZOD 
ZOE 
2079 


Seuisr " 


"Have you created & stared your oor 
Ms eekitsaPp Cy Gales Mie eo STM set 
IF At="N" THEN 1126 
IF AB="Y" THEN 1140 
SIT 14a 
Peer © DIS © LISP " Please refer ta the" @ OLSr 


Bi=e 
Biome 
Genes 7a 
DiIsP @ DISP "Enter problem suboroaqram 
UN Sl i ae 
Eeeen@e Disk @€ ani Sr 
Poo SRAPHIIC 
eee eo ise " a) 
IF N=2 THEN FLAG 11 
RETURN 
SFLAG 
CALL 


MOAT USER MANUAL" @ 
iene s Gi Mie Gl am Sac. 


a 


file name. 


select 
DISELAY" 
LATA UISRPLAY" 


pragress caption..." @ OLSP 


@ INFUT N@ CFLAG 11 


eee it Si 
Wer ALT ” 


Sl ocmes vate tay UUM ae 


Geral Set seem net: “MEET” € CALL “ACU 


Le 
TALL 


Cuenca a ot | 
"GRAD" @ CALL 


laadinga PGi. ae . 


SRR Ome oii ile inh Opie iS 
eS 
Eee 
ITFLAG 
RE TIIRN 


Satrse “Almast Tifisned.«.< 
"NEWH" @ CALL "RONV" @ CALL 
9 @ CLEAR 


Pe Riis? @ (CALL. “TERM: 


na) aa) 





SMOY) 
SOLA 
SAL 
SEV 


(LEAR @ DISP USING 
IMAGE 2/7,12X%,48 
Oia USING aose » Ni,N2 


3 
IMAGE 1/,4X, 20," ign 


2010 5 Pt 


lies Variables",/,4X%,20," Inequality 


“ins craints” 


3040 
=) = (2) 
HAD 
3074 
tt 
OEY) 
2A'7 iD 
s199 
sila 
SL2o 
3120 
SQa0) 


tse 


4010 
Ad29 
4On 


to the 


4040) 
4050) 
4aE 
4070 
CY 
4029 
4079 
410%) 
4119 
$120 
41320 
414M 
Sy 
ra Bait") 
4140 


ta be changed" 


417 
4152 
4170 
4200 
42710 
A229 
4229 
7 or 
4240 

Zug 
4240 


ap Nz; ENTER)" 


DIGF USING toeS@ +. Fo 

IMAGE 2/,1X,"“Initial Llesian: F = 
IF N2=0 THEN zaze 

IF AMAX(G)<C2(3) THEN Dt="FEAQSIEBLE” 


ae 


Peso e = ThNrPens Leese 


Wiag tletNG “L1x, 1A” = fe 
IN KREY# 1, "EDIT" GOSUB 4aa@ 
ON KREY# 4, "CONTINUE" GOTO 2130 


mer LAEEL 

ae Sle 

RETURN 

Peeee © LISse © WIAP “Select editing 
fete eaish YORTArLeS" 

eee se ols 2s ChUERAUL TT VAaElEesa' 

MIN N GOTO 4@30,4140 

CLEAR @ DISP "Enter the address 
changed" @ INFUT I 

Peo ISIN 4050 ° 1, Xacl) 

Pimeise 2/,1%X, "current ia tds a ee Fs 

Mere o lice “Enter the new value..." @ 

Pea=on © LISP € nise ag bata eae) Goo 7" 

@ INFUT At 

THEN 4000 

THEN 4114 


REC Were), eo wll ae 


INPUT 


Lohan mo ote Variable 


it ) = iM : by 
INFLIT X@(T) 
2 01SF " 


IF At="N" 
ice ai= i" Y it 


GAT 4070 

1, Wee 

Wem + 6) 2 KOC) FY, dG Lt), OD 

Gate 230 

Peesneecettoe CRIS “select default array..." @ DIS &€ 
. Wee h Invegers ” 

eae eC niSsr " eee) pea | Sy eS ENE LIT Ty 

CLEAR @ DISP "Enter the address in C of the variable 


@ INPUT I 

IF N=1 THEN N@=C1iCI) ELISE NoaClecl) 

IF N=1 THEN Gt="Cic" ELSE E¢="C2c" 

BisrP USING 4200 +: ES,1,N® 

DMewe 2/7,1k, "current value: ",3A,K,°) = 
NISP @ DISP "Enter the new value..." @ 
IF N=1 THEN CicI)=N® ELSE CSCI) =Na 
EEEGR @ OFISP @€ DISF “Editing camp] 
N; ENTER)" @ INPUT At 

IF At="N" THEN 4000 

Ir At="Y¥" THEN 424 

RE TURN 


ak 
INFLIT N@ 


ete @ © DIS ” ( 


=i 





Module: LOGO 


Called by: Autost 


Bugectaonm: Generate MDOT "welcome" graphic display on CRT. 


Seratched from memory upon execution. 
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1 She "Loca" 

1@ PEN 1 @ GCLEAR 

20 SCALE @,10,0, 14 

3@ FEN -1 @ GCLEAR 7 

40 PEN 1 @ GCLEAR = 

5@ CSIZE 2® @ LORG 5S @ FEN -1 

40 MOVE 5,5 @ LABEL "mi" 

7@ (SIZE & @ MOVE 5.4.5 @ LABEL "x" 
30 WAIT 1900 @ FEN 1 @ GCLEAR 

99 CSIZE 12 @ LORG = 

19@ MOVE 1,2 @ LABEL "m" 

110 MOVE 1,4 @ LABEL "qd" 

12@ MOVE 1,4 @ LABEL "a" 

120 MOVE 1,2 @ LABEL "t" 

14@ CSIZE 2 
15@ MCIVE 2,2 
140 MOVE 2,4 
17@ MCIVE 2,4 
120 MOVE 2,2 
170 SUBEND 


iD () fe fa 


Pebet Serwenmrlter pase d ” 
LABEL “esign” 

Pease fT tinleac tam 

YS ES aac a Mie 


Do DH 


aS 





Module: DEFALT 


Called by: Autost 


Function: Set MDOT working parameters to their default val- 


ues. DEFALT is called each time a new problem is 


executed. 
Nomenclature: 
SG) Maximum number of iterations 
en G2) Consecutive iterations for convergence criteria 
C1(3) Variable metric method select: 0 = DFP, 1 = BFGS 


C1(4)-C2(10) Unassigned 


C2 (1) Finite difference perturbation factor 

e2(2) Minimum absolute finite difference step 

e203) Violated constraint criterion (tolerance) 
C2(4) NeGilve constrains eriterion (thickness) 

25) Push-off factor multiplier (theta zero) 

e200) Maximum value of push-off factor 

C2(7) Factor used in DIRECT when infeasible 

C2(8) Factor used in step length estimate based on F 
C2(9) Factor used in step length estimate based on X 
6210) F convergence criterion (relative) 

C201 1) F convergence criterion (absolute) 


C2( 12) Defined zero 
S213) Epsilon, used to prevent division by zero 


C2(14)-C2(20) Unassigned 
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oO 


wer 





E2(2 


1 @vY) 


Sle hee ply” 

METICN BASE 1 

Pte treo SF INTEGER NI,N2,Na,N4,01,W2,03,04,C1018) 
@) 

Be ue l=Zer 1) 

MAT CB=ZER( 20) 

Mic¢ids2O ' max # iteratians 

meee y=2 ! onset. Canv. 


feces y—l §' win H update 


Peete. © Fin. Diff. mult. 
Pee = Hi “min, Tin. ster 
C2(S)=.904 § canst. viol. 
meee v=-. i ft active canst. 


Sees =1°! pish—aff mutt. 


fees § max. push-att 
Pengo 17000 § Pri (OIRECT) 

Beees=. ! aby. mult. 

Pees y=. ! des. var. mult. 
RmSC1@)=.001 $ min rel F 
Peet y=.80) § min abs F 
Beets )=.081 |! sere 

fs =.8OO1 ! enasiian 
SLIREND 


Spal 


5 ote F 





Module: 


Called by: 


Function: 


Module: 


Called by: 


Function: 


PROB 


User 


Provide the user with a skeleton subprogram into 
Which the optimization problem is entered by 


editing. 


PROB (edited) 


Autost, GRAD, FDSRCH, VMSRCH 


Input the number of design variables and the 
number of inequality constraints (90-100). 
fopuGeene Intial design and side constraints 
(130-210). 

Evaluate the objective function and inequality 


constraints (230-9000). 


Nomenclature: 


K 1 
L$ 
US$ 
K1-X0 


Flag to indicate first call or subsequent call 
String used in lower bound input 
String used in upper bound input 


Design variable names used in problem input 
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i@ &s 


Meme Chi, %0),F,¢).L¢).010)) 


SeemeTIuN EASE 1 


A 


44 [i 


=a) 
7a F 


Lae 
1a5 
130 
140 


(72) 
PAQD 
7010 
7FaLO 
‘FNAO 
‘70S 2) 
WAH 
FOED 
7070 
POS 
FO'70 
7100 


im. (= 


Spesieimd AL, AL, XS,X4,X5,X6,X%7, X38, X7, 


CM PSC4) , INTEGER N1,N2,N2,N4, 01,02, 003, 04,0101) 


1) 


IM L#C41,U8C4) 
F K1>1 THEN 22 
EAD N1i,N2 

DATA 


Gy 


feet —2 THEN SUBEXIT 


FOR =i TO Nit 
READ X(1) 

IF N2=@ THEN 1! 
READ Lt, it? 


fees = NN THEN LCI) 


TF Uet="N" THEN 
NEXT I 

TIATA 

DATA 

CATA 

DATA 


= DATA 


LATA 
DATS 


: DATA 


LIATA 
DATA 
HSE yapta 


7% 


Sate eee ee VALA LS) 
Web vate ey? Stee ui Pye ula.) 


' liser-defined expressians 
Beeoeeecive Tanction 


WE=Sla+ f 

f CONSTRAINTS 
SUBENT 
XL=X¢1) 
XS=SX(2) 
Aa=X CS) 
X4=X(4) 
X2=X (5) 
XO=X04) 
X7=X(7) 
AS=X CE) 
X'7=XC7) 
XM=X(190) @ RE 


PM MM HM i. M 9 fp 
= 
“Tl 


N1=7 
TURN 


THEN 


= THEN 
ES alan 


THEN 
Vinleis 


THEN 


THEM 
THEN 
THEN 


RETURN 
RETURN 
RE TIIRN 
RE TUR 
RE TIIRN 
RETURN 
RE TURN 
RETIIRN 
RETURN 


see 





Module: COON TY 


Called by: Autost 


Calls: TERM, PROG, ACON, GRAD, DIRECT, FDSRCH, CONV 


Funetion: Control constrained optimization by the method of 


feasible directions (Fig. 2). 
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1O SUB "CCUNT" (X@0),F9,G0),L0),u0)) 

20 OPTICN BASE 1 

30 COM F$C4] , INTEGER N1,N2,N2,N4,01,02,02,04,C1(10) , SHORT 
C2(20) 

4@ SHORT X1(19),X(18),0010),G1(12,10),AB(12),5¢10) 

SQ SHORT F,E0,AQ,F1 

G® INTEGER A1(12) 

70 IF FLAG() THEN Tape 

@ REDIM X1i(Ni),00N1),G1(N1+2,N1),S¢(N1) 

JQ A®,F1,E@=0 @ F=FO 

100 MAT S=ZER(N1)@ MAT X=x@ 

110 MAT AL=ZER(N1+2)@ MAT Az=A1 

12® Gi=Mi+1 

130 IF GisCi(1) THEN 140 

140 SFLAG 1 @ CALL "TERM" ( X(),F,G(0) ) @ SUBEXIT 

15 SUBEXIT 

160 CALL "PROG" ( X@(),X(),F/F@,F ) @ IF FLAG('7) THEN ZUBEXI 


Pee ae AUN Cat), ALC),Ae() ) @ IF FLAG(Y) THEN SUBEXIT 
Poin Cerne ALE VTERM" ( X();F,G(0) ) 
Boge Fenis( ig) TREN SLIBEXIT 


240 TF Na+N4#0 THEN REDTIM GitCN3+N4,N1) 

eee Nei), amit) tee, Abt ,Ase(),Git,) ) @ IF FLAG 
Peet MeN SUBEXIT 

220 IF Nia+N4#® THEN 240 

220 MAT S=-0@ GOTO 240 

Seca tlie) Cl CEC) ECD, AL(), ASI), S(),B80 ) @ IF FLAG(Y 
) THEN SULBEXIT 

foo ErPLAG S 

2460 IF MAXABCS)SC2(12) AND ABS(EQ)tCe(12) THEN 240 

279 3FLAG 3 

Soe Net THEN CALL "TERM” ( X(),F,0() ) 

Soir Fee (ide THEN SUGEXIT 

Semmiieeigewnme che) =Ge(l2) THEN CALL "TERM" ( X(),F,G¢) ) 

eae ir FLAB(Ii®) THEN SUBEXIT 

Bee ee ee lei) hUTHEN GALL UT TERM" € X(),F,G0) 3} 

Seer PeAais(i@) OR FLAG(T) THEN SIUBEXIT 

240 C2(4)=C2(4)/3 

Sarit 170 

260 MAT S=C1L/MAXAB(CS))F#5@ CFLAG 4 

370 MAT X1=X@ FI=F 

ee eee emeotenwiim € X(),LC) uC) FSO), &C) AL), ASL), GIC,), S¢ 
),AY ) 

SO Fiepist?7) THEN SUBEXIT 

404 IF FLAG(4) THEN 289 

410 CFLAG 3 

429 MAT X=(1)FX1+(AQO) FS 

eee Cer, AC), F ,uC) ) 

ame emmemimNe a X(),X1(),F,Fi > © 
430 IF FLAG(S) THEN CALL "TERM" ( X 
eee) itR FLAGS) THEN SUE 
470 ‘SLUBENT 


IF FLAG('7) THEN SUBEXIT 
(igs oad) ) 
EXIT ELSE GOTO 120 


SS 





Module: UCONT 


Called by: Autost 


Calls: GRAD, PROG, TERM, VMSRCH, CONV, NEWH 


[muMecton. “control unconstrained optimization by the 


variable metric method (Fig. 1). 


Nomenclature: 
K Indicates an iteration in which no move parameter 
to improve the design was found 


A One-dimensional search move parameter 
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10 SUR "CONT" (xO) ,FO,G0),L0),0I0)) 

20 OPTION BASE 1 

SQ COM F$C6] , INTEGER N1,NZ,N3,N4,01,02,02,04,c1(ia) 
C2 (26) 

4Q@ SHORT X(10),X1(1@), 0010), 01(19),3¢10),H(1%, 10) 

S® SHORT F,F1,A8 

6@ INTEGER E 

7@ IF FLAG(9) THEN SUBEXIT 

SO RECIM X(N1),X1¢N1),000N1), D1 (N21) 

70 MAT X=X9@ F=F9 @ K=a 

19@ CALL "GRAD" ¢ Xt),F,0),G¢) ) 

11@ MAT H=IDN(N1,N1) 

12@ G@i=t1+1 

13@ CALL "PROG" ¢ x@(),X0),F/FO,F ) @ IF FLAG(7) THEN SUBEXI 


HIRT 


meee tr Ui-ig1(1) THEN iva 

PeQeeceimiat © HALL “TERM” ¢ X(),F,50) ) 

Peele Plast id) Ge FLAGCY) THEN SUBEXIT 

1798 MAT S=HFE 

120 MAT s=-% 

170 MAT Z=(C(1/MAXAEB CS) ) #5 

<9 MAT X1=X@ MAT DL=O@ FLi=Fr 

21@ CALL "YMSRCH" ¢ X€),F,00),50),8 ) @ IF FLAG(?) THEN SUBE 

may 

meer Are etis) THEN 278 

230 KS +] 

S@0ekr hie THEN 11 

Peeeeeieis + 8 TALL "TERM" ¢ X¢),F,6¢) ) 

Seomimereping 1) WR FLAGYY) THEN SLIBEXIT 

279 MAT X=(1)#X+(48) F535 

Meee eee CO XX) SAL) CF, 1l ) @e IF FLAG(Y) THEN SUBEATT 

fever FlAi(sa) THEN CALL "TERM" ¢ eo ae ) 

320 TF FLAGCI@) OR FLAG(Y) THEN tl 

eegte FPLARCG) THEN 119 

ee eee maa < XC), FP, BC) iO) ) © IF FLAGSY) THEN SUBEXIT 

emer CXC) KIC) CC) I), y 2 © TF FLAGCY) THEN 
Boece xX i] 

S40 miiyTo pee 

350 SLUBENT 


Si 





Module: 
: Called by: 


Function: 


ACON 


CCONT 


Determine the number of currently violated (N3), 
and active (N4), inequality constraints. 
Construct the Ajl() vector of the addresses in G() 
of the active/violated constraint set. 
Consuructethe Az2() vector of the current values 
Ciel s. Se t. 

Violated constraint information is stored in the 
first N3 rows of A1l() and A2(), active constraint 


information in the last NH rows. 
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ese TRIN (is€),A1(),AzZC)) 
Boon tIN EASE 1 


Seieemer el Ga) [INTEGER Mi,N2@,Na,N4, 01,082, fra, @4,C1(1o) 


It (28) 

40 CIN KEY# 1 GOT sao 

ar PLAis(7) THEN SUBEXIT 
£0 IEFLAG 2 @ Na, N4=a 

momeam f=) Fu No 

Seige (s(f) eS s) THEN 118 
70 IF GCI) 2=n2(4) THEN N4=N4+1 
190 GOTO L2H 

mo Ns=N St] 

heer NEAT I 

120 [TF NSS =Ni*\2 THEN Lae 
Psa SPLAG 2 2 SLIBEXIT 

120 ITF NS+N4=0 THEN SLIBEXIT 
140 MAT AL=ZERCNSB4+N4)@C MAT AZ=AI 
179 Ws, 14=1 

meee [=] Tr Nz 

fotete Intl) e=Celt(s) THEN 230 
By TP G1 )ACe (4) THEN 230 
“19 ALCNS+14) =I 

S20 AD(NS+4)=GC1) 

2 WI4=1g4+] 

=49 GOTO 290 

ma el (a p= 

260 Ax CHA) =O0T) 

ey tas lot] 

Bec NEXT I 

290 GOTO 210 

200 ‘FLAG Y @ 'BUBEATT 

S10 SUBEND 
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Module: GRAD 


Called by: CCONT, UCONT 


Gals: PROB (edited) 


Function: Calculate the gradient of the objective function 
Dyetinsce forward finite difference approximation. 
In constrained optimization, calculate gradients 


of the active/violated inequality constraint set. 


Nomenclature: 


F2 Intermediate function evaluation 


N Design variable perturbation 


10 








19 SUE "GRAD" ¢X0).F,00),60),010),AZ0),G1(,)) 


Seep pig FASE 1 


30 COM ECS] , INTEGER N1,N2,NE,N4, 01,02, 02,04,01(19) 


C2 (20) 
40 SHORT F2,X7(1@),N 
S® ON KEY# 1 GOTO 210 
é£@ IF FLAG(?) THEN SUBEXIT 
76 FOR I=1 TO NI 
SO MAT X7=X 
7O N=C2(1)*ABS(X9(1)) 
100 IF NECE(2) THEN Neca 2) 
110 X7CI)=X7C1)4+N 
foe CALL PP ( 2.%94),F 2,60) 
126 IF NZ=@ THEN 1280 
140 IF NZ+N4=0 THEN 1&0 
SQ FOR J=1 TO NZ+N4 


160 Gitd, TD SCa CAL CI) KAS CUIO/N 


170 NEXT 

ower lL) =(F 2=r)/N 

BVO NEAT I 

200 GTO 22a 

eee srlAl ¥ S&S SUEBEAIT 
220 SUBENE 


ZN 


4 


= 


' 
t 


aS 


7 
4 





Module: 


Called by: 


DIRECT 


CCONT 


Bumetron: oOlvye the direction-finding subproblem in the 
method of feasible directions. 
Calculate constraint push-off factors (170-200). 
Initialize working arrays for currently feasible 
(2000-2080), or infeasible (1000-1090) designs. 
Determine the direction vector S(), and the para- 
meter BO (300-700). 
Nomenclature: 

ae, Working array constructed from G1(,) and T() 

B(,) Working array initialized as -A’A 

BO Kuhn-Tucker parameter 

Be Intermediate element value used in pivoting 

B3 Intermediate element value used in pivoting 

B9 Intermediate variable used in pivoting 

ce) Working vector 

DOC) Working vector initialized as D() 

GO(,) Working array initialized as G1(,) 

G9() Working vector used in constructing G0O(,) 

T9() Working vector of element indices 

J9 Working integer scaler 

K9 Working integer scaler 

N9 Dimension of active/violated constraint set 


me 





Et) Working vector 


EO) Vector of constraint push-off factors 
UC) Working vector 
Y() Solution vector, contains S() and BO 
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19 SUB "DIRECT" (04),G1(,),A10),020), 50), Be) 
20 OPTION BASE 1 

3@ COM FECG] , INTEGER Ni,N2,N2,N4,01,02,05,04,01(19) , SHORT 
2 (20) 

4@ SHORT B2,B2,89,A(12,11),7¢12),B(13,13),P¢11),U613), C013), 
Y¥(11),0@(1) ,GO(12,10),G7(10) 

S® INTEGER J9,K97,17(12),N9 

40 ON KEY# 1 GOTO 450 

79 IF FLAG(Y) THEN SUEBEXIT 

30 NI=N3+N4 

90 REDIM T(NY),PCNi+1), ¥(N1+1) 
1a@ FOR I=1 Ta Ny 

11@ T(1)=Ce(S)*(1-A2(1) /02(4) 72 
120 IF T(1)2C2(4&) THEN T(1)=C2(4) 
130 NEXT I 

140 MAT FP=ZER(Ni+1)@ MAT D@=D 

15@ MAT [i@=(1/MAXAE (DO) )#ra 

140 MAT Gd=G1 

17@ FOR I=1 Ta NY 

180 MAT GI=G0(1,) 

199 MAT G9=(1/MAXAB (GY) )#GY 

200 MAT GO(I,)=G7 

210 NEXT I 

220 IF NS?® THEN GOSUB 1900 ELSE GOSUE zaaa 
230 MAT EB=A#TRN(A) 

240 MAT B=-E 

250 N=SUENDCB, 1) 

240 MAT IY=ZER(N) 

270 AI, KI=0 

230 FOR I=1 Toa N 

290 IF C(I) >=@ THEN 23a 

300 BY=C(1)/BC(I,1) 

210 IF B9<=A7 THEN 230 

220 K7=1 @ AT=EP 

420 NEXT 1 

349 IF K=O THEN 540 

350 I9=17(K2) 

B40 L9(KI) =O 

370 IF J9=0 THEN ID( EF) SKI 

320 BOER(EZ, KT) 

39% FOR I=1 To N 

40@ B(KY,1)=B(K9,1)/B2 

419 NEXT I 
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$20 CCK) =AP 

430 ECE, E7)=S=1/B2 

440 FOR T=1 TO N 

330 IF [=K? THEN S20 

4490 BS=BCI,ET) 

470 BCI,ERY)=a 

430 FOR J=1 TO N 

Peemet lt J =EC LT) I) -Batetk7, J) 
=90 NEXT J | 

s10 M¢T)=0 CI) -B3#A9 
geeroMEXT I 

23M ISNT 27Y 

=40 FOR I=1 TaN 

359 lICI)=a 

zoo ey =17¢( 1) 

a7 OofF J72.=o THEN s'7o 
eo UC TIaC(Iy) 

270 NEXT I 

£40 MAT Y=TRN(A) Fl 

419 MAT Y=P-Y 

620 MAT S=¥C1:N1) 

630 BO=YC(N1i+1) 

640 GATA 440 

ea srlLAb 7 @ SUBEXIT 
£60 SIIBEND 

1aoHW MAT C=ZER(NY)@ MAT L=ZERCN'Y) 
1919 MAT A=CIN(NY,N1I+1) 
1929 MAT AC, 1:N1)=G8 
1930 MAT AC,Nit+1i1)=T 

1949 MAT FCIIN1L)=Do 

1420 MAT F=-P 

19640 FCNI+1)=2¢7) 

1970 MAT l=AF 

19s0 MAT tC=-C 

1970 RETURN 

“000 MAT C=ZER(N7+1)2 MAT US=ZERCNT+1) 
2010 MAT A=KON(N7+1,Nit+1) 
2929 MAT ACLENY, 1EN1)=00 
2930 MAT ACLINY,Nit1l)=T 
“9044 MAT ACNY+1,1:N1)=00 
2030 FCNI+1)=1 

2068 MAT C=A(,N1+1) 

24970 MAT C=-£ 

=980 RETURN 


vg= 





Module: 


Called by: 


Calls: 


Function: 


FDSRCH 


CCONT 


PROB (edited) 


Perform one-dimensional search for constrained 
Sperm zarvionm im the method of feasible 
directions. 

Estimate initial search move parameter (1000- 
1180). 

Check for side constraint violations (2000-2050). 
Establish bounds on solution, feasible (3000- 
3340) or infeasible (4000-4360). 

Retime  solucion DY polynomial approximation, 


feasible (5000-5240) or infeasible (6000-6580). 


Nomenclature: 


Move parameter 

Working vector of move parameters 

Working vector of move parameters 

Initial A based on change in objective function 
Initial A based on attaining feasibility 

Array of constraint values during search 
Porencoauotnot DO) and) S@) or Git{i,) and S() 
Working vector of constraint gradients 

Working vector of maximum constraint values 


Working vector of objective function values 
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Poeeteeeiered” (€X¢).L(),110),F,0(),60),81(),820),614,),80),A 


20 OPTICN BASE 1 

30 IMOM F204] , INTEGER N1,N2Z,N3,N4,01,@2,02,04,C01(19) , SHORT 
Be (2) 

$9 SHORT X7°10),R1,R2,A1,A2,AC4), 764) ,404), 804, 30) ,00,0G0010) 
,PO,P1,FP2,BO,Aa(S1) 

SONI EGER I,K 

40 ON FEY? 1 GOTO 210 

Por FLAGi 7) THEN SUBEXIT 

S0 REDIM X7(N1),B¢4,N2),G@(Ni), A@(NZ+1) 

Pueeie=( soA (a) )/S © R2=s-Rl 

19@ A,AL,AZ,AS,AC1)=40 2 CFLAG 4 

110 MAT M=ZER(4)@ MAT X7=X@ Y(1)=F @ MAT &C(1,)=52 MC 1) =AMARX ( 
a) 

120 GUSLIE Loan 

120 ACBI=H=A @ GOEWE Boome 

1490 ITF Na=4 THEN GOSUE seGge EL 
150 [IF N3s=0 THEN GOSUE See EL 
140 A=ACE) 

meee Acie (i2) THEN lo 

130 SFLAG 4 @ SUBEXIT 

1970 FHeYcK) @ MAT GHEECE, ) 

200 GOTO 229 

Pees “7 C SUBEXIT 

fo SIRE MEA 

1900 GOSHE Lease 

1010 IF Nst#0 THEN lasa 

1920 A=81 @ RETURN 

1930 IF Azrz*AL THEN A=2*A1 ELSE A=MAX (CAL, AZ) 
190496 RETURN 

1930 DM=00IT (0, 3) 

19440 AL=C203)F AES (CF) /ABS (DO) 

1070 FOR f=1 Ta Nil 

IWEO A=C2(7)FABSCXCI)) /SABSCSCI)) 

1a70 IF Atal THEN AL=A 

mov NEXT I 

1119 IF NS=0 THEN RETURN 

PiegePun T=1 Th Ns 

1120 MAT G@=G1¢I,) 

1144 DO=ONIT (8, >) 

1150 A=-(AZ(I)/DO) 

114@ IF A*SAZ THEN AZ=A 

1174 NEXT I 

liso RETURN 

2800 MAT X=(C1L)#FX9+(A) F5 

2419 FOR I=1 To Ni 

Seco mie AC L)AL¢(L) THEN X¢(T)J=t¢1) 

203@ IF X¢(IT) UCI) THEN X¢(T)=UICT) 

2940 NEXT I 

249540 RETURN 


=. 
} 
ae’ b 
~~ 
= 
—_ 


FRC 4¢ae 
E GOENSHIE Aaa 


Uo 





“= @) 
2010 
O20 
AQ) 
24) 
SAE 
3960 
2070 
OV) 
3a 
A100 
211) 
2120 
SLIM 
314a 
21S 
2160 
el 
21S 
s1590 
2Z2On 
3710 


3220 (2) 


bike & 


MAT BC(2, 


M(1)=AMAX (5) 
Meuse C2 KA) 
Mi) =AMAXK (13) 
YC) =F 


JF,G0) ) 


)=0 
Pomeece ds -1 C1) THEN 3070 
Te MCs) =f2ati2) THEN S31! 
ACS)=HACLIIFRIF(ACE ")-A(1)) 
A=A(3) 
mOSuUB TOaM 
ele be « 2X) 
en) 
eee) = 
MAT Gt) =G 
RETURN 
ne 
ee = yk 
a) 
MAT G=B(2, ) 
MAT B(3,)=6 
AC(ZI=H(1+R2) FACS) -REFA(1) 
A=A(2) 
GOSUB 2aan 
elie Foe 
Y(2)=F 
M(t) =AMAX (5) 
mer + Ete , 8 
ee Ce et () 
die Mi) SCZ 12) 
AC1)=AC) 
sec t= (3) 
MCL) =MCS) 
MAT G=E(2, ) 
MAT B(1,)=6 
GOTO Sis 


=) 


JF,G0) ) 


eS COee. sor. 


TREN SRE Tinh 
THEN RETURN 


Viz 





42) 
4010 
Ai) 
BARE) 
Aad 
mS bal) 
406% 
4a7a 
4020 
4070 
4100 
4119 
4120 
41°30 
4140 
4139 
$140 
4179 
419@ 
41°70 
4204) 
4210 
4220 
42238) 
4240 
4250 
426.0 
4279 
Peis) 
4270 
40) 
4310 
22 
49 
424@ 
47350 
4360 


M(1)=AMAX (1G) 

meme Pt € 2,X0),F, G0) 2 
¥(2)=F 

M (2) =AMAX (13) 

MAT E(2, )=06 

IF M(2)>M(1) THEN 4080 
IF M(2)2@ THEN 4140 

IF ¥C2)<¥(1) THEN 4140 
ACB)=A(1)+R1#(A(2Z)-AC1)) 
A=A(3) 
GOS Baa” 
CALL PS ( 
¥(S) =F 

M (2) =AMAX (13) 

MAT BCS, )=G 

RETURN 

AZ) =A(2) 

YCR)=¥ (2) 

MAT G=B(2,) 

MAT BC3,)=6 

MCS =M(2 

AC 2)=(14+R2) #A (3) -RE*AC1) 
A=A(2) 
ISHShle Soa 
CALL FR ( 
Y(2)=F 

M( 2) =AMAX (3) 

MAT E(2,)=6 

IF M(2)2M(C3) THEN RETURN 
IF M(2)20 THEN 4310 

IF ¥(2)2¥(3) THEN RETURN 
A(1)=A(3) 

WCL)=EYCR) 

M(1)=MC 2) 

MAT G=E(3,) 

MAT B(1,)=6 

GOTG 4140 


i ae a. 


pee er cere i Oe ae | 


739 





SMEG) 


pee eGo lo, (ACs SACL I I= CYC 2)-Y¥(1))/ CAC2)-A01579/708 


Regent.) ) 


2410 
Sa) 
SN @ 


aH @) 


PL=CY(CBI-Y¥ CLI) /S (ACE) -ACL) )-FEH#(QC1)4002)) 
AM(1L)=-(P1/ (ZEPB)) 
FOR I=1 TO Nz 


P2=((BC3,1)-B(1,1))/( ACB) -AC1) )-(BC2, ID-B(1, 1) / (A024 


mo) / AC S)—AC2)) 


ATA 
Ai) 
=a 7 @ 
a] OO Pett 2) 
aAT iy) 
a1ag 
—a119 
as Ok) 
el 2) 
2140 
ie ci) 
—aLéw 
Sly vt 
Sis 
=190 
mt AY) 
21% 
me) 
St ow) 


ta 4) 


PI=(B(2,1)-B(1,1))/(AC2Z)-AC1) PRE (ACL) 4A(2)) 
PO=B(1,1)-P1#AC1)-P2#A(1)°2 

EQ=P 1" 2-4#POsP2 

IF B@>® THEN 5110 

A@(I1+1)=-(PO/FL) 

GOTO Size 

AOC I+1) =MAX (-F1+S0R (BO), -PL-SOR (EO) ) 

AOC T+1)=AO(I+1)/(2#F2) 

NEXT I 

A(4)=AMIN(A®) 

A=A(4) 

BOSE [Oa 

CALL FS ( 2,X0),F,G0) ) 

Y(4)=F @ MAT E(4,)=5 

M(4)=AMAX (G) 

k= 1 

FOR I=2 Ta 4 

IF ¥(1)<¥CK) AND 
NEXT I 

RETURN 


Hii =620be) THEN K=1 


Sie) 





6000 MAT AQM=LER(CN+1) 

40035 IF AMIN(M) +0 THEN 4:24@ 

paeme— 0 C3.) —-¥(1))/( ACS) -AC1) Y-CY¥ (2)-¥ 01))9/7(802)-8(1)))708 
eA) ) 

G020 FLISCYCSZI-¥C1LII/ CACSZI-AC1) D—FP2F(8Q¢01)4002) ) 
4420 AMCLT)J=—-CPI/ (SRF) ) 

6ag@ FOR T=1 TO Na 

4050 GOs 449% 

é060 TF Beso THEN 4a72 

£O70 AOCT+1)=-(FO/F 1) 

442” GOTO 4110 

4070 ABOaOCT+1)=MIN(-FPitsoR (BO), -PLl-SORCE@) ) 

4149 AOC T+41)=AG( IT +1) / 02°F 2) 

eoig NEXT I 

4120 ACI) =AMAXK (CA®) 

6120 MAT A@=ZER(CNZ-N2)@ J1f=1 
S140 FOR T=1 To N2 

4150 FOR J=1 Tl No 

4140 IF T=Ai(.t) THEN 4240 
6170 NEXT wt 

41330 GOSUEB 4540 

6170 IF Base THEN 464220 

4200 AOC IL)=-(CFO/FI1) 

4210 Ji=J1+1 @ GOTO 4248 
4220 AO(ILI=SMIN(-Pi+S0R (EO) ,-Pi-sQR CEQ) ) 
6220 AOCJLIHAGCIL)/(C2RF 2) 
42:33 J1=Ji¢+i 

i240 NEXT I 


=n 





* a2 
AAO 
S279 
apne 2) 
coe 7 U) 
LOLA) 
431% 
AA 


ACS) =MIN(CAC4),AMIN(C AD) ) 


A=A(4) 

KOSUE Soon 

CALL FS ( 2,X(),F,G() ) 
Y(4)=F @ MAT B(4,)=6 

M4) =AMAX (15) 

i= 

FOR I=2 To 4 


350 IF ¥Y¥CIJAYCR) AND MCI)<=02(04) THEN k=] 

&t4QO NEXT I 

4320 RETURN 

eee eee Ca SCl) 7 ACS) ACT JK I(MC2)-—M6 1) / (ACS) -AC1)))/ 64 
Se) >A( 2) ) 

ee? OF I= (MC a1) 7 CAC 2) AL 1) )-P2tCAC1)+AC2)) 


(28) 
6370 
44 0@ 
aoe 
6420 
£420 
6440 
E450 
£448 
470 
44:30 
&4'7@) 


AC4)=-(P1L/ (SRF) ) 

IF ACS)<A(4) AND AC4)2 ACS) THEN 6420 
M4) =29M(2) 

GUT 64460 

A=A (4) 

ISCrslle Saad 


Print 25% (),F, 150) ) 

M(4)=AMAX(G) @ ¥(4)=F @ MAT E(4,)=6 

M=AMIN(M) 

K=AMINRGW 

RETURN 
Pz=((B(S,A1(1))-B(1,A1(1)))/ (ACS) -AC1)) -(B (2, 4161) )-E Ct 


sALCIDIISCAC2ZKACL) DDS CAS) -AC2)) 


G50 
ae, 
6352) 
ABO 
a=48) 


Owe) / (ACs) -A(Z 


& S58) 
4548 
457% 
robe bath 2) 


P1=(B(2,A1(1))-EB(1,A161)))/(ACZ)-AC 1) )-PE* (ACL) +A(2)) 

PO=B(1,A1(1))-P1#A(1)-P2#A( 1) 2 

RO=P1“2-4#PO*P2 

RETURN 

P2=((BCR,I)-B(1,1))/(AC3)-A(1) = (BCE, 1)-B(1, 1) )/(ACZ)-A 
Dy) 

Pis(B(2,1)-B(1,1))/(ACZ)-8(1) )-P2*( ACL) +002) ) 

PO=B(1,1)-P1#A(1)-PE¥A(1) 72 

EQ=P1°2-4#PORP2 

RETURN 


Sa 





Module: 


Called by: 


Calls: 


Function: 


VMSRCH 


UCONT 


PROB (edited) 


Perform unconstrained one-dimensional search in 
variable metric method. 

Esmabiasne boumds on the Solution (100-350). 
Refine the solution by polynomial approximation 


(360-520). 


Nomenclature: 


AO 
AC) 
BC) 
DO 
FQ 
ae) 
rer.) 


Move parameter 

Working vector of move parameters 

Working vector of move parameters 

Dot product of D() and S() 

Working value of objective function 

Working vector of objective function values 


Working vector of objective function values 


8 3 





19 SUE "VMSRCH" (X(),F,00),5¢),0@) 
20 OPTION BASE 1 

30 COM FSC4] , INTEGER N1,N2,NZ,N4,01,02,032,04,01(19) , SHORT 
C220) 

4@ SHORT [I0,A1(11),X9(10),¥(4),A(4),R1L,R2,P1,P2, 2,164), F9, 
B 

S@ ON KEY# 1 GOTO S4a 

40 IF FLAG(?) THEN SUBEXIT 

7@ RI=(B-SOR(S))/2 @ RB=(1+S0R(5))/2 
20 MAT A=ZER(4)@ MAT Y= 

70 VC1)=F @ F9=0 

190 D@=poaT (Oo, Ss) 

11@ MAT A1=ZER(N1+1) 

120 MAT X'3=X 

1320 AL(1)=C02(S)#ABS(F) /ABS (D0) 

140 FOR I=1 TO Ni 

150 AL(I+1)=C209) #OBS(X9(1))/ABRCR(I)) 
14% NEXT I 

17@ A(Z)=AMIN(AL) 

120 MAT X7=(1)8#X9+(A(2) ) #8 

170 CALL P# ( 2,X9(),F9 ) 

200 VYCR)=F9 

210 IF Y(2)>Y¥(1) THEN 210 

230 ACB)SA(2) @ YC3)=¥(27) 

230 AC 2)=(1+RZ)#A(2)-R2#A(1) 

240 MAT X3=X 

250 MAT X9=(1)#X9+( AZ) ) #5 

240 CALL PR ¢ 2,X70),F7 ) 

270 YC2)=F7 

220 IF Y¥CR)2YC3) THEN 244 

29@ ACL)=ACS) @ ¥(1)=¥(3) 

200 GOTO 22a 

S1@ ACB)=R1I*A(2) 

320 MAT X7=X 

3SO MAT X9=(1)#X9+(A(3) ) #5 


a4 





340 CALL PR ( 2,X70),F9 ) 

S50 ¥(S)=F9 

SOQ PR=(YCR)-¥C1))#(A(2)-AC1)) / (ACB) “ACL HON C2) HV 01) 4 0ACR)- 
ACi))) 

37@ PR=PR/ (ACE)-AC1) +(ACR)-AC2) )400)/((A(2)-A(1)) #0A03)-AC1) 
y#(AC3)-ACZ))) 

SQ@ PZ=C (OY CZ)-¥C1))/(ACZ)-AC1) )-0@) / (ACE) -AC1) )-PS#(244(1) +8 
ea) 

290 PLslid-2eP2#A( 1) -S#PRHAC 1) O2 

400 R=P22-2#P1#PS 

41@ IF E>=0 THEN 420 

420 J=z @ GOTO 470 

430 AC 4)=(-F2+S0R(B))/(S#PS) 

440 MAT X7=X 

ASO MAT XP=(1)#X7+(A(4)) #3 

440 CALL P# ( 2,X90),F9 ) 

A7® ¥(4)=F9 

420 J=4 

49@ MAT Y1=¥C 12) 

Sao F=AMIN(Y1) 

S10 K=AMINROW 

S20 AQ=A(K) 

S20 GOTC She 

540 SFLAG 7 @ SUBEXIT 

S50 SUBENT 


eS 





Module: NEWH 


Called by: UCONT 


munmetion. Update the approximation to the inverse of the 
Hessian matrix used in determining the search 


direction in the variable metric method. 


Nomenclature: 


DOC ;) Update matrix 

EO) Working vector initialized as X()-X1() 
5 Dowmproduct of FC) and Y() 

ri) Working vector 

T10) Working vector 


eC .,.) Working array 
te, Working array 


YO) Working vector initialized as D()-D1() 


86 





meme "NEWH" (X<(),X1(),0¢),01(),H!,)) 

20 OF TION BASE 1 

0 COM F#C4] , INTEGER N1,N2,N2,N4,01,6:2,03,04,C1¢19) , SHORT 
C2(20) 

4 SHORT 019(1@,10),¥(1),P (10), 7T1(19),7T2614,19),7T(1),8, 72419 
,10) 

SO CN EEYH 1 GOT ze 

40 IF FLAGC7) THEN SUBEXIT 

7a REDIM D@(N1,N1),/(N1),P¢(N1),TLCN1),T2ZCN1,N1),TSCN1,M1) 
a0 MAT F=X-X1 

7@ MAT Y=DI-[t1 

100 S=O0T(P,Y) 

110 MAT T1i=H#Y 

120 MAT TZ=TI*TRN(T1) 

120 MAT T=TRN(Y)#T1 

140 MAT DO=FP#TRN(P) 

150 MAT D@=( (54+T(1) #0103) ) /S°2) #00 

149 MAT DO=(1) 80+ ((C103)-1)/T01) 0 #T2 

170 MAT TZ=TI#TRN(F) 

180 MAT TRB=F#TRN(T1) 

170 MAT T2=T2+TS 

20@ MAT [@=(1)#0@+¢-(01¢C3)/5))4T2 

210 MAT H=H+re 

220 GOTO 24@ 

230 SFLAG 7 @ SUBEXIT 

244 SUBEND 


=) 





Module: CONV 


Called by: CCONT, UCONT 


Function: Determine whether the design has converged to the 
optimum in the last iteration. 
Update convergence criteria based on iteration 


hes POR y. 


88 





1@ SUB "CONV" (X(),X10),F,F1) 

20 OF TICN BASE 1 

3Q@ COM PEC4] , INTEGER N1,N2,N3,N4,01,02,03,04,01(19) , SHORT 
2 (20) 

40 SHORT X9(10) 

SQ IF FLAG(Y) THEN SUBEXIT 

40 CFLAG 4 

7@ C2(3)=(02 (38) +0BS((F1-F)/F1L))/2 

#@ MAT X9=X1-X@ MAT X7=X9/X1 

YO C207) 5 (C209) +MAXAB( XY) )/2 

100 IF ABS(F1-F)<MIN(C2(13),C2(11)*ABS(F1)) THEN @3=034+1 ELS 
E Ga=@ 
119 IF ABR(F1-F)/MAX(ABS(F1),.@@@@1)<C2(10) THEN G4=04+1 ELS 
E W4=0 

120 IF MAX(G2,04)>=C1(2) THEN SFLAG & 
120 IF MAX(02,0:4)>@ THEN SFLAG 4 

14@ SUBEND 


Ss 





Module: PROG 


Called by: CCONT, UCONT 


Bunmecron: Generate optimization progress information 


OUcCDUL. 


90 





SUE "PROG" (X@0),X0),Y,P) 
OPTION BASE 1 
COM FS04] , INTEGER N1,N2Z,N2,N4,01,02,03,04,C1¢19) 


2) ) 


UN REYSH 1, "INTRET" Gora 3200 


Hoe eee.) THEN SUBEXIT 
IF FLAG(11) THEN 230 
IF G@171 THEN 1388 


BECERR © FEN 1 @ LORG & 


SCALE -2, 21, “3.2, 3, 2 
XAXTS @,1,9, 20 
YAXI'S d, ies ote 
FUR i=—3 To 3 
iaives—.o,i © LAEEL I 
NEXT I 
MIBVE 16,3 

Mive 10), 2 


Boceee iteration Mistir y” 
me? 
MEV =a, 1, 


AGieis s 


@ 
7 
See Poet be ler Ty Cerrunt 


FOR I=1 TO Nt 

MOVE G1,X(I)/X@(1) @ LAFEL I 
NEXT I 

MOVE C1,Y @ LABEL "fF" 
SUBEXIT 

IF G@1+1 THEN 274 


CLEAR @ DISP WSING "12xX%,40" : Ft 
Bese oi veratian Kh jective Function" 
Per LAEEL 

eer als VN be 7 EEN S 


IMAGE 2X, 20,12X,K 
OLITE Bile 

SFLAG 9 @ SUBEXIT 
SLUBEND 


=) J. 





Module: TERM 


Called by: CCONT, UCONT 


Pemetion; Generate output of optimization results. 


92 





ieee Temi: (XC) ,F,&C)) 

“4 HPTIOIN EASE 1 

34 COM FP'Cé) , INTEGER N1,N2,N3,N4,01,02,03,04,C1(019) , SHORT 
ee (2) ) 

4@ IF FLAGCY) THEN SUBEXIT 

20 Inve 1aad 

60 IF FLAG(1) THEN GOS 2aae 

70 TF FLAG(C2) THEN GOSUB saae@ 

= TF FLAG(CS) THEN GOSuUB 4aaea 

70 IF FLAG(4) THEN Gosue Saae 

194M TF FLAG(S) THEN GOSUIE Amaa® 

1190 IF FLAG(11) THEN 15 

120 PRINT @ FRINT @ PRINT 

1:34 GRAPH @ TOFY @ SUBEXIT 

l4o FRINT © FRINT @ Cory 

120 ‘SUGENLI 

140” FRINT @ FRINT 

Pee PRING WelNis "12k,48" = FS 

mOce PRING USING tao: F 
meso MAine 1/,2%, OPTIMUM 
i“#4@% FOR T=1 Ta Ni 
PeeomPREN, USING twee + 1,X¢(1) 

moeemeyoe 17,10K," KC" Ky") = "LE 

107% NEXT I 

i@s® IF Ne=@ THEN 11230 

197o FUR T=1 To Ne 

meormeNt USING 1ii@ ; 1,G¢1) 

meromrniee 17, 1ax, "isc", hk," ) = “Sk 

Pic aE XT if 

meow aN)! © PRINT @ PRINT "Termination tased an:" @ PRINT 
1140 SFLAG 14 

115@ RETURN 

BVO rmRiNt URkcee¢ded nan. mo. of iteratians.” @ RETURN 

3400 PRINT “Excessive number of violated mince tnosa  C 

RETURN 

4000 FPRINT "Failure ta find a direction te i imprave the desi 
gn." @ RETURN 

S000 FRINT "Failure ta find a mave parameter ta imprave the 
design." @ RETURN 

6000 PRINT "“Canvergence"” @ RETURN 


ud ; kz 


oS 





Listings of the edited versions of PROB which were used 
for the test cases presented in Chapter IV are included as 
examples of MDOT problem input. The subprogram used to enter 
the unconstrained test problem was renamed "BANANA", while 


that edited for the constrained problem was renamed "BEAM", 


94 





Pars 
2a 0 
10 i: 


pg Oi 


UB "EANANA" 
FTIUN EASE 1 
eee & J 
(4) ) 


Gi ( ) lcm) 


, INTEGER N1,N2,N3,N4,01, 02,03, 04,119) 


SS SHORT X1,X2,X2,X4,X5, X6, X7, XB, X9, XO 


49 0D 
aig I 
7@ R 
19@ 
1@> 
138 
148 
142 


IM L#04],U8043 
F K1>1 THEN 220 
FAD N1,N2 

DATA 2,0 


IF F1=@ THEN SUBEXIT 
FOR I=1 TO Ni 
Rem XC] ) 


IF N2Z=0 THEN 17 


a Peay el Ba 


Bi =Va Gli 


Vee Gr — lB ese ICL p=VAL (lle) 


meee AL SPOR KL EEK SETS KS 24h 1 2a 2eXith 


RETURN 
RETLIRN 
RETURN 
RE TILIRN 
RETURN 
RETURN 
RETURN 
RETWRN 
RETURN 


144 READ L#, it 

130 IF L?="N" THEN L(I)= 
140 IF Ue="N" 

170 NEXT I 

eet ATA -1 

£92 DATA 1.5 

<9 DATA 

=94 DATA 

2035 DATA 

29& DATA 

297 DATA 

ZO: DATA 

297 TATA 

=~190 DATA 

22% ISQSUE 7O1a 

2290 !' User-defined expressians 
fom a yective Function 
4a) 

477 W2=Ws+l 

“00 ! CONSTRAINTS 

740% SUBENT 

7@1M@ Xi=X(1) @ IF Ni=1 THEN 
FHLO XS=X(2) @ IF Ni=2 THEN 
YO2Q X3=XC3) @ IF Nl=sa THEN 
2M@4W X4=X(4) @ IF Ni=4 THEN 
YOaO@ X3=X(5) @ IF Nl=S THEN 
TOED X&=X(4&) @ IF N1l=4 THEN 
FQ70™ X7=X(7) @ IF Ni=7 THEN 
F7OSO XS=X(S) @ IF Ni=s THEN 
7OIO XV=X(7) @ IF N1l=? THEN 
7100 KO=X(19) @ RETURN 


95 


SHORT 





19 SUE "BEAM" (K1,X0),F,G0),L0),010)) 

20 OF TION BASE 1 

30 COM PSC6] , INTEGER N1,N2,NS,N4,01,0:2,03,04,C1(19) , SHORT 
2 (20) 

SS SHORT X1,X2, XS, X4, XS, X4,X7, XS, KF, XO 


49 DIM L#CS1,Ue(4] 
S@ IF Ki>1 THEN 220 
79 READ Ni,Nz 


19% DATA 2,3 
1935 IF KLl=0 THEN SUGBEXIT 
ime tit (Ni) ,INCNI), ACN1),GCNZ) 


FOR T=1 TO Ni 
REAL X¢I) 

IF N2=0 THEN 170 
READ Lt, Us 
eee ON 
TIF Uet="N" 
NEXT I 
DATA 3. m 
DATA 14, 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
ATA 
GOSUB 2PoHLa 


130 
149 
178 
ZO4 
OZ 
203 
O04 
ZO 
ZOE 
297 
SAS 
207 
210 


220 


anti wtie =e ELSE (1) SVU ) 
TREN CI =) see ect ce AIC 1 SVrt Cole) 


oY 4 
1,20 


230 ! User-defined expressians 

3201 G=X1l 

242 H=X2 

B77) Wejyestive function 

400 F=ZOOrkheH 

477 O2S2+1 

29a !' CONSTRAINTS 

S01 G¢1)=400/ (B*H*2)-1 

=92 6(2)=106444.7/ (BFH™S) -1 

202 GC2)=H/14-B 

79ow SUBEND 

7@19 Xi1=X(1) @ IF Ni=i1 THEN RETURN 
7920 X2=X(2) @ IF Ni=2 THEN RETURN 
7920 X3=X(3) @ IF Ni=3 THEN RETURN 
JO40 X4=X(4) @ IF Nil=4 THEN RETURN 
7O@590 XS=X(5) @ IF Ni=S THEN RETURN 
7040 XS=X(S&) @ IF N1l=4 THEN RETURN 
907@ X7=X(7) @ IF Ni=7 THEN RETURN 
7WOsO@ XS=X(S) @ IF Ni=a THEN RETURN 
JOIO X7=X(7) @ IF Ni=7 THEN RETURN 
FiDQ XG=X(10@) @ RETURN 


Se 





10. 


tee 
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